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GLOBAL  FREQUENCY  DISTRIBUTION  OF 
EXOSPHERIC  TEMPERATURE 


Jack  W.  Slowey 


INTRODUCTION 

It  was  desired  to  determine  the  distribution  of  exospheric  temperature 
as  given  by  certain  models  of  the  earth's  upper  atmosphere.  The  required 
distributions  were  to  cover  all  conditions  of  solar  and  geomagnetic  activity 
as  well  as  spatial  location  above  the  earth.  They  were  to  be  employed  to 
facilitate  the  use  of  the  models  in  general -perturbation  orbit  integration. 

The  models  of  interest  were  the  1971  Jacchia  model  (Jacchia,  1971)  and  the 
model  given  for  heights  above  120  km  in  the  U.S.  Standard  Atmosphere  Supple¬ 
ments,  1966  (COESA,  1966).  For  convenience,  these  models  will  be  referred 
to  as  J71  and  US66,  respecti vely . 

The  problem  was  approached  in  three  stages.  In  the  first,  the  historical 
record  of  solar  activity  over  the  past  12  solar  cycles  was  used  to  establish 
the  effect  of  solar  activity  on  the  desired  distributions.  In  the  second,  the 
effect  of  the  diurnal  temperature  variations  as  given  by  the  models  was  imposed. 
Finally,  the  distribution  of  geomagnetic  disturbance  was  estimated,  and  its 
effect  on  the  exospheric  temperatures,  as  given  by  the  models,  was  taken  into 
account.  These  three  stages  of  the  investigation  and  the  results  from  each 
are  described  in  the  following  sections. 

SOLAR-ACTIVITY  VARIATION 

The  earth's  upper  atmosphere  undergoes  a  large  variation  in  temperature 
and  density  in  the  course  of  the  11-year  cycle  of  solar  activity.  The  range 
of  this  variation  is  not,  however,  constant  from  one  cycle  to  the  next  but 


varies  in  accordance  with  the  level  of  solar  activity  at  maximum,  a  highly 
variable  quantity.  Thus,  it  is  necessary  to  take  a  large  portion  of  the 
historical  record  of  solar  activity  into  account  in  order  to  determine  the 
statistical  effect  on  the  distribution  of  atmospheric  temperature. 

The  historical  record  of  sunspot  numbers  goes  back  to  1749.  However, 
in  most  models  of  the  thermosphere  and  exosphere,  including  those  under  con¬ 
sideration  here,  the  variations  associated  with  solar  activity  are  correlated 
with  the  10.7-cm  radio  flux  from  the  sun.  The  record  of  the  10.7-cm  solar 
flux  goes  back  only  to  1947  and  is  inadequate  for  the  present  purpose.  It  is 
necessary,  therefore,  to  be  able  to  relate  the  observed  sunspot  numbers  to 
the  corresponding  10.7-cm  flux  for  use  in  the  atmospheric  models.  The  rela¬ 
tion  that  was  used  was 

Fio  7  =  49.4  +  0.97  R  +  17.6  exp  (-0.035  R)  ,  (1) 

where  7  is  the  mean  adjusted  10.7-cm  solar  flux  and  R  is  the  mean  Zurich 
sunspot  number.  This  is  a  very  slightly  modified  version  of  an  equation 
previously  developed  by  L.  G.  Jacchia  and  the  author.  It  is  reported  (Euler, 
Lundquist  and  Vaughan,  1978)  to  give  a  correlation  coefficient  of  0.98  with 
a  data  base  of  F^g  7  from  67  to  260,  where  the  means  of  F-jg  7  and  of  R  were 
13-month  running  means  of  the  monthly  mean  values  (the  13-month  running  mean 
is  centered  on  the  month  in  question  and  gives  half  weight  to  the  months 
6  months  before  and  6  months  after  that  month). 

The  distribution  of  the  nighttime  minimum  in  the  exospheric  temperature 
was  obtained  by  sampling  every  third  monthly  value  of  the  13-month  running 
means  of  the  sunspot  number.  Smoothing  over  6-month  intervals  is  usually 
recommended  in  connection  with  the  atmospheric  models.  The  13-month  means 
which  are  routinely  used  in  solar  activity  prediction  -  were  already  available, 
however,  and  it  was  clear  that  the  slightly  greater  smoothing  would  have 
very  little  effect  on  the  resulting  distribution.  Equation  (1)  was  used  to 
obtain  the  mean  10.7-cm  flux,  which  was  then  converted  to  temperature  by 
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Table  1.  Distribution  of  global  minimum  exospheric  temperature  as  computed 
from  the  J71  and  US66  atmospheric  models. 


Temp.  Interval  (°K) 

J  7 1  Model 

US66  Model 

n/nt(  ) 

f(K_1) 

N 

n/nt(x) 

fU"1) 

0 

- 

+  10 

64 

12.0 

1.20  x  lo'2 

61 

11.5 

1.15  x  io"2 

+  10 

- 

+  20 

45 

8.5 

0.85 

41 

7.7 

0.77 

+  20 

- 

+  30 

34 

6.4 

0.64 

32 

6.0 

0.60 

+  30 

- 

+  40 

24 

4.5 

0.45 

23 

4.3 

0.43 

+  40 

- 

+  60 

34 

6.4 

0.32 

33 

6.2 

0.31 

+  60 

- 

O 

CO 

4- 

27 

5.1 

0.25 

29 

5.5 

0.27 

+  80 

- 

+  100 

39 

7.3 

0.37 

32 

6.0 

0.30 

+  100 

- 

+  120 

36 

6.8 

0.34 

33 

6.2 

0.31 

+  120 

- 

+  140 

32 

6.0 

0.30 

29 

5.5 

0.27 

+  140 

- 

+  160 

29 

5.5 

0.27 

24 

4.5 

0.23 

+  160 

- 

+  180 

22 

4.1 

0.21 

30 

5.6 

0.28 

+  180 

- 

+200 

23 

4.3 

0.22 

19 

3.6 

0.18 

+200 

- 

+230 

23 

4.3 

0.14 

26 

4.9 

0.16 

+230 

- 

+260 

19 

3.6 

0.12 

27 

5.1 

0.17 

+260 

- 

+300 

22 

4.1 

0.10 

19 

3.6 

0.09 

+  300 

- 

+350 

23 

4.3 

0.09 

26 

4.9 

0.10 

+  350 

- 

+400 

13 

2.4 

0.05 

15 

2.8 

0.06 

+400 

- 

+450 

8 

1.5 

0.03 

12 

2.3 

0.045 

+450 

- 

+550 

12 

2.3 

0.023 

13 

2.4 

0.024 

+550 

- 

+650 

3 

0.6 

0.006 

6 

1.1 

0.011 

+650 

- 

+  750 

0 

0.0 

0.000 

2 

0.4 

0.004 

Relative  to  the  temporature  corresponding  to  R  =  0  (Fin  7  =  67).  This 
temperature  is  596.1  K  for  the  J77  model  and  603.2  K  for  the  IJS66  model. 

+N  is  the  number  of  data  points  in  the  interval,  NT  is  the  total  number  of 
data  points  (=  532). 
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J 7 1  MODEL,  T0  =  379°0  +  3°24  F 


4 


Figure  la.  Distribution  of  global  minimum  exospheric  temperature  over  12+  solar  cycles, 
J71  model . 


US66  MODEL 


Figure  lb.  Distribution  of  global  minimum  exospheric  temperature  over  12+  solar  cycles, 
US66  model . 


using  the  appropriate  equation  in  the  particular  model.  Sampling  began  with 
the  value  for  December  1847  and  ended  with  that  for  September  1980,  giving  a 
total  of  532  data  points.  The  11  points  at  the  end  were  taken  from  the 
National  Oceanic  and  Atmospheric  Administration's  predictions  of  July  1978. 

The  resulting  distributions  for  the  two  models  are  tabulated  in  Table  1  and 
are  shown  plotted  for  the  J71  and  US66  models  in  Figures  la  and  lb,  respectively. 
The  relative  frequency,  f,  is  the  fraction  of  the  total  per  1-degree  temperature 
interval.  The  model  equation  used  to  relate  the  nighttime  minimum  temperature 
to  the  mean  10.7-cm  flux  is  given  on  each  of  the  figures.  The  mean  values 
resulting  from  these  distributions  are  726  K  for  the  J71  model  and  748  K  for 
the  US66  model.  The  mean  10.7-cm  flux  corresponding  to  the  mean  exospheric 
temperature  is  approximately  107  for  both  models. 

DIURNAL  VARIATION 

The  diurnal  variation  in  exospheric  temperature  as  given  by  either  the 
J71  or  the  US66  model  can  be  written  as 

T^TqU  +  RD)  ,  (2) 

where  To>  is  the  exospheric  temperature,  TQ  is  the  nighttime  minimum  exosphere 
temperature,  1  +  R  is  the  (constant)  ratio  of  the  maximum  to  the  minimum  global 
exospheric  temperature,  and  D  is  the  diurnal  parameter.  The  diurnal  parameter 
varies  between  0  and  1  and  is  given  by 

D  =  sinm  0  +  (cosm  n  -  sin"1  o)  cos11  (j)  , 

with 

o  =  2  U  +  4>g  I 

o  -  2*  I  'l*  ~  'f’g  I  » 

T  =  H  +  3  +  p  sin  (H  +  y)  ,  (3) 
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where  <f>  is  the  latitude  of  the  point  in  question,  is  the  latitude  of  the 

maximum  in  exospheric  temperature,  and  H  is  the  hour  angle  of  the  sun.  In 
the  case  of  the  J71  model,  <|>R  is  set  equal  to  the  declination  of  the  sun, 
while  in  the  US66  model,  ij>g  =  0.  Constants  for  the  two  models  are 

R  =  0.3  ,  m  =  2.2  ,  n  =  3.0  ,  3  =  -37°  ,  p  =  6°  ,  y  =  43° 

and 

R  =  0.28  ,  m  =  1.5  ,  n  =  2.5  ,  3  =  -45°  ,  p  =  12°  ,  y  =  45° 

for  the  J71  and  US66  models,  respectively. 

A  computer  program  was  written  to  integrate  numerically  the  area  on  the 
globe  between  the  point  of  minimum  exospheric  temperature  and  the  exospheric 
isotherm  defined  by  a  particular  value  of  the  diurnal  parameter.  Results 
from  this  program  are  tabulated  for  both  models  in  Table  2  and  are  plotted 
in  Figure  2a.  Table  2  gives  results  for  the  two  extreme  cases  of  equinox  and 
solstice  in  the  case  of  the  J71  model.  As  can  be  seen,  the  differences 
between  the  results  for  the  two  cases  are  quite  small.  Thus,  the  results  for 
the  single  case  of  equinox  were  taken  to  represent  the  distribution  for  any 
time  of  the  year  in  the  subsequent  computations  with  the  J71  model. 

The  diurnal  distribution  of  exospheric  temperature  is  quite  different  in 
the  two  atmospheric  models  as  a  result  of  the  differences  in  the  model  constants. 
This  is  evident  from  the  isotherm  plots  included  with  the  published  models. 

It  can  also  be  seen  in  Figure  2b,  where  the  diurnal  temperature  distribution  — 
which  depends  on  the  slope  of  the  diurnal  area  curve  in  Figure  2a  —  is  plotted 
for  each  of  the  two  models  for  the  case  where  the  minimum  exospheric  temperatures 
are  equal  to  the  mean  values  given  above.  In  the  J71  model,  the  frequency  is 
much  larger  near  minimum  temperature  than  it  is  at  higher  temperature  because 
the  isotherms  are  much  more  widely  spaced  near  minimum  temperature  than  they 
are  for  higher  temperatures.  In  the  US66  model,  the  variation  in  spacing  is 
much  less,  and  the  isotherms  are  actually  closer  together  at  minimum  than  they 
are  at  maximum. 


Table  2.  Integrated  global  isotherm  areas  as  a  function  of  the 
diurnal  temperature  parameter,  D,  for  the  J71  and  US66 
atmospheric  models  (total  area  =  1). 


D 

J71  Model 

♦B  =  0?0  *B  =  2  3°  5 

US66  Model 
*B  =  020 

0.00 

0.0000 

0.0000 

0.0000 

0.05 

0.1037 

0.1014 

0.0397 

0.10 

0. 1787 

0. 1748 

0.0841 

0.15 

0.2452 

0.2398 

0.1307 

0.20 

0.3063 

0.2995 

0.1791 

0.25 

0.3635 

0.3551 

0.2287 

0.30 

0.4174 

0.4074 

0.2794 

0.35 

0.4685 

0.4576 

0.3310 

0.40 

0.5158 

0.5061 

0.3831 

0.45 

0.5625 

0.5533 

0.4357 

0.50 

0.6055 

0.5991 

0.4886 

0.55 

0.6473 

0.6437 

0.5414 

0.60 

0.6883 

0.6870 

0.5938 

0.65 

0.7286 

0.7290 

0.6454 

0.70 

0 . 7685 

0.7694 

0.6965 

0.75 

0.8078 

0.8090 

0.7475 

0.80 

0.8469 

0.8481 

0.7983 

0.85 

0.8855 

0.8866 

0 . 8490 

0.90 

0.9239 

0.9248 

0.8996 

0.95 

0.9620 

0.9626 

0.9500 

1.00 

1.0000 

1.0000 

1.0000 
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URNAL  PARAMETER, 


MODEL,  Tft =  726 


Global  distribution  of  exospheric  temperature  from  J71  and  US66  models  computed  for  mean 
value  of  global  minimum  exospheric  temperature. 


The  global  distribution  of  exospheric  temperature  was  calculated  for  the 
two  models  from  the  corresponding  distribution  of  minimum  temperature  by 
computing 


fG<Tl'  V 


-£ 


WAVA2  ~  AP 


T2  '  T1 


(4) 


where  fg(T-|  ,  T^)  is  the  global  relative  frequency  per  degree  in  the  interval 
between  and  T2,  fg^)  re^atlve  frequency  per  degree  of  the  minimum 

temperature  Tg,  and  and  are  the  global  isotherm  areas  corresponding  to 
the  two  diurnal  parameters  computed  from 


D. 

1 


(5) 


The  sum  was  taken  with  steps  of  1°  in  T^  and  values  of  A  were  computed  by 
interpolation  from  Table  2. 

The  resulting  distributions  (K  =  0)  are  given  for  both  models  in  Table  3 
and  are  plotted  in  Figure  3a  for  the  J71  model  and  in  Figure  3b  for  the  US66 
model.  Mean  values  of  the  exospheric  temperatures  given  by  these  distributions 
are  817  and  855°K,  respectively,  for  the  J71  and  US66  models. 

GEOMAGNETIC  VARIATION 


To  determine  the  effect  of  the  geomagnetic  variation  on  the  distributions, 
it  was  first  necessary  to  obtain  the  distribution  of  geomagnetic  disturbance. 
The  historical  record  of  the  geomagnetic  index  commences  in  1932  and,  hence, 
encompasses  the  last  4  complete  solar  cycles  (17-20).  Only  data  beginning 
midway  through  cycle  19  were  readily  available  to  the  author,  however.  Data 
from  cycle  19  were  excluded  both  because  they  represented  only  a  partial  cycle 
and  might  be  biased  for  that  reason  and  because  cycle  19  itself  was  extremely 
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Table  3.  Global  distribution  of  exospheric  temperature  from  the  J71 
and  US66  atmospheric  models  for  K  =  0. 


J71  Model  US66  Model 


★ 

Temp.  Interval  (° 

K) 

n/nt(*) 

f(°K_1) 

n/nt(%) 

f ( °k" 1 ) 

0  - 

+20 

2.8 

0.141  x  io"2 

1.2 

0.061  x  10 

+20  - 

+40 

4.6 

0.229 

2.7 

0.134 

+40  - 

+60 

4.9 

0.247 

3.6 

0.180 

+60  - 

+80 

5. 1 

0.254 

4.3 

0.217 

+80  - 

+  100 

5.5 

0.274 

5.0 

0.251 

+  100  - 

+  120 

6.0 

0.298 

5.6 

0.282 

+  120  - 

+  140 

6.3 

0.315 

6.2 

0.309 

+  140  - 

+  160 

6.6 

0.328 

6.7 

0.333 

+  160  - 

+  180 

6.7 

0.335 

6.8 

0.340 

+  180  - 

+200 

6.1 

0.303 

5.6 

0.281 

+200  - 

+230 

7.2 

0.241 

7.1 

0.236 

+230  - 

+260 

6.2 

0.205 

6.5 

0.215 

+260  - 

+  300 

7.1 

0.178 

7.6 

0.191 

+300  - 

+  350 

7.1 

0.141 

7.8 

0.155 

+350  - 

+400 

5.2 

0.104 

6.1 

0.122 

+400  - 

+450 

3.8 

0.075 

4.6 

0.092 

+450  - 

+500 

2.8 

0.055 

3.5 

0.070 

+500  - 

+550 

2.1 

0.041 

2.6 

0.052 

+550  - 

+600 

1.5 

0.030 

2.0 

0.040 

+600  - 

+700 

1.6 

0.016 

2.4 

0.024 

+700  - 

+800 

0.7 

0.007 

1.2 

0.012 

+800  - 

+900 

0.3 

0.003 

0.6 

0.006 

+900  - 

+  1000 

0.1 

0.001 

0.2 

0.002 

+1000  - 

+  1100 

0.0 

0.000 

0.1 

0.001 

+1100  - 

+  1200 

0.0 

0.000 

0.0 

0.000 

«r 

Relative  to 

596?  1 

K  for 

•  the  J71  model  and  603?2  K  for 

the  US66  i 

model . 

Table  4. 

Distribution  of  the  K( 
recent  solar  cycle.  11 

geomagnetic  index  over 

the  most 

K 

P 

atg*  (°K) 

N 

f  =  n/nt 

0 

0.0 

180 

.0448 

.0448 

0+ 

9.3 

254 

.0632 

.1080 

1- 

18.8 

286 

.0712 

.1792 

1 

28.1 

309 

.0769 

.2561 

1+ 

37.4 

315 

.0784 

.3345 

2- 

46.9 

351 

.0874 

.4219 

2 

56.2 

333 

.0829 

.5047 

2*- 

65.6 

304 

.0757 

.5804 

3- 

75.2 

316 

.0786 

.6590 

3 

84.6 

333 

.0829 

.7419 

3+ 

94.1 

260 

.0647 

.8066 

4- 

103.9 

207 

.0515 

.8581 

4 

113.5 

157 

.0391 

.8972 

4  + 

123.5 

127 

.0316 

.9238 

5- 

134.0 

93 

.0231 

.9520 

5 

144.4 

61 

.0152 

.9671 

5+ 

155.4 

41 

.0102 

.9774 

6- 

167.5 

29 

.0072 

.9846 

6 

180.1 

23 

.0057 

.9903 

6+ 

194.1 

10 

.0025 

.9928 

7- 

210.4 

9 

.0022 

.9950 

7 

228.9 

5 

.0012 

.9963 

7+ 

251.0 

3 

.0007 

.9970 

8- 

279.0 

4 

.0010 

.9980 

8 

313.4 

4 

.0010 

.9990 

8+ 

357.6 

2 

.0005 

.9995 

9- 

417.5 

1 

.0002 

.9998 

9 

495.1 

1 

.0002 

1.0000 

X  -  28“ 

K  +  0?03  exp 
P 

(Kp). 
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Table  5.  Global  frequency  distribution  of  exospheric  temperature  from  the 
J71  and  US66  atmospheric  models  with  the  geomagnetic  temperature 
variation  included. 


J71 

Model 

US66 

Model 

Temp.  Interval  (  K) 

n/nt(x) 

F( °K~ 1 ) 

N  ( X ) 

F(°K-1) 

0  - 

+20 

0.2 

0.010  x  io~2 

0.1 

0.004  x  10'2 

+  20  - 

+40 

0.8 

0.040 

0.4 

0.020 

+40  - 

+60 

1.6 

0.082 

0.9 

0.047 

+60  - 

+80 

2.5 

0.127 

1.6 

0.081 

+80  - 

+  100 

3.4 

0.171 

2.4 

0. 121 

+  100  - 

+  120 

4.2 

0.211 

3.3 

0.163 

+  120  - 

+  140 

4.9 

0.243 

4.0 

0.202 

+  140  - 

+  160 

5.4 

0.268 

4.8 

0.238 

+  160  - 

+  180 

5.8 

0.288 

5.4 

0.269 

+  180  - 

+200 

6.0 

0.301 

5.8 

0.289 

+200  - 

+230 

9.1 

0.302 

8.8 

0.293 

+230  - 

+260 

8.5 

0.282 

8.3 

0.278 

+260  - 

+  300 

9.8 

0.244 

9.8 

0.245 

+  300  - 

+  350 

9.8 

0.196 

10.2 

0.204 

+350  - 

+400 

7.6 

0.152 

8.3 

0.166 

+400  - 

+450 

5.8 

0.116 

6.6 

0.132 

+450  - 

+  500 

4.3 

0.085 

5.1 

0.102 

+500  - 

+550 

3.2 

0.063 

3.9 

0.078 

+550  - 

+600 

2.3 

0.046 

3.0 

0.059 

+600  - 

+700 

2.8 

0.028 

3.7 

0.037 

+  700  - 

+800 

1.3 

0.013 

2.0 

0.020 

+800  - 

+900 

0.6 

0.006 

1.0 

0.010 

+900  - 

+  1000 

0.2 

0.002 

0.5 

0.005 

+1000  - 

+1100 

0.1 

0.001 

0.2 

0.002 

+1100  - 

+  1200 

0.0 

0.000 

0.1 

0.001 

Relative  to  596 ? 1  K  for  the  J71  model  and  603?2  K  for  the  US66  model. 
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US66  ATMOSPHERIC  MODEL 


eluding  geomagnetic  variation,  US66  model 


m 


unusual.  On  the  other  hand,  the  data  from  cycle  20,  which  coincides  very 
nearly  with  the  mean  of  cycles  8-20,  seemed  idealy  suited  to  determination  of 
the  distribution.  This  was  done  by  sampling  the  Kp  index  once  each  day 
over  the  interval  of  cycle  20.  The  resulting  distribution  is  given  in 
Table  4  and  plotted  in  Figure  4.  Figure  4  also  shows  the  distribution  of 
the  increase  in  exospheric  temperature  that  results  from  the  derived  Kp  dis¬ 
tribution.  The  temperature  increase  is  the  same  in  both  models  and  is  given 
by 


ATg  =  28°Kp  +  0?03  exp  (K  ) 

The  distribution  of  geomagnetic  temperature  increase  was  applied  to  the 
global  distributions  for  K  -  0  in  a  manner  similar  to  that  used  to  apply  the 
diurnal  variation.  The  relative  frequency  corresponding  to  a  given  tempera¬ 
ture  increase  was  determined  by  interpolation  of  the  values  given  in  Table  4. 
The  resulting  global  distributions  are  given  for  both  models  in  Table  5  and 
are  plotted  in  Figure  5a  for  the  J71  model  and  in  Figure  5b  for  the  US66 
model.  Mean  values  of  the  exospheric  temperature  from  these  distributions 
are  881  K  for  the  J71  model  and  918  K  for  the  US66  model. 

CONCLUSION 

The  frequency  distribution  of  exospheric  temperature  taken  over  all 
locations  on  the  earth  and  all  conditions  of  solar  and  geomagnetic  activity 
was  derived  for  both  the  J71  and  US66  atmospheric  models.  Both  distributions 
have  a  fairly  sharp  maximum  at  a  temperature  just  above  800  K.  The  skewness 
of  the  distributions  is  reflected  in  somewhat  higher  values  of  the  weighted 
mean  temperature.  The  mean  temperature  is  approximately  880  K  in  the  case 
of  the  J71  model  and  920  K  in  the  case  of  the  US66  model. 


Subsequently,  distributions  were  determined  using  all  8  3-hourly  values  of 
the  K  index,  both  with  and  without  the  cycle  19  data.  These  were  considerably 
smoother  near  the  maximum  in  the  distribution  but  not  significantly  different 
in  any  other  way. 
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